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A spin-i magnetic impurity coupled to a one-dimensional correlated electron system have been 
studied by applying the density renormalization group method. The Kondo temperature is sub- 
stantially enhanced by strong repulsive interactions in the chain, but changes non-monotonically in 
the case of electron attraction. The magnetization of the impurity at zero-temperature shows local 
Fermi liquid behavior. 
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During the last decade, much effort has been devoted 
to the investigation of the electron-correlation effects. 
They are present most prominently in the copper-oxide 
materials which display high-T c superconducting and un- 
usual normal-state properties and may even show 
heavy-fermion behavior In one dimension (ID), an 
interacting electron system is known to be a Luttinger 
liquid J3|. It seems possible that in the future by means 
of quantum wire fabricated by nanotechnique one can 
study the effect of magnetic impurities on Luttinger liq- 
uidg. 

The problem of a Kondo impurity coupled to a Lut- 
tinger liquid was first considered by Lee and Toner ||. 
By using perturbative renormalization group theory they 
found that the Kondo temperature Tk depends on the 
coupling J between the impurity and the conduction 
electrons in the form of a power law (not exponential) 
provided that the interaction between electrons is suf- 
ficiently strong compared to J. Because of the neglect 
of local backscattering their treatment does not preserve 
SU(2) symmetry. The local backscattering was included 
within the poor man's scaling by Furusaki and Nagaosa 
H who obtained a stable strong coupling fixed point for 
J > as well as J < 0. Moreover, by employing a 1/J 
expansion, they found that a temperature expansion of 
thermodynamic quantities shows critical behavior, i.e., 
non-integer exponents. It is an open question though, 
how well such a perturbative expansion is justified in the 
strong coupling regime. In fact, a simplified model 
for conduction electrons, which consists of right-moving 
spin-up and left-moving spin-down electrons, leads to the 
conclusion that the Kondo impurity behaves like a local 
Fermi liquid . On the other hand, boundary conformal 
field theory methods ||[l0) result in two types of critical 
behavior, i.e., either a local Fermi liquid or one of the 
form proposed in ||. Most recently, Chen et al Jll| pre- 
dicted local Fermi liquid behavior by making use of the 
parity and spin-rotation symmetry of the problem. 

The above literature indicates that the critical behav- 
ior of a Kondo impurity coupled to a Luttinger liquid is 
not fully understood yet. One scenario would be that of a 
local Fermi liquid but with a substantially enhanced Tk 



due to strong correlations between conduction electrons. 
For example, the characteristic energy scale in the heavy 
fermion system Nd2-xCe K Cu04 is much larger than ex- 
pected for Kondo ions coupled to free electrons fl2|| . 

In this paper, we study ground-state properties of a 
Kondo ion antiferromagnetically coupled to a ID Hub- 
bard model by applying the Density Matrix Renormal- 
ization Group (DMRG) Q. The Hamiltonian is: 



H n = - 



cLcj 



-1,<T 



h.c] + U ni\hii + JS ■ sq 



where S is the impurity spin, Sq an electron spin at i = 
and other notations are standard. In the limit of U = oo 
and half-filling, the problem reduces to one of an impu- 
rity coupled to a Heisenberg antifcrromagnetic chain, a 
system studied before (see Refs. Jl4]-^6|). Here we will 
consider all U values and band fillings equal to as well as 
different from one-half. 

The Kondo problem is concerned with the question of 
how the magnetic moment of the impurity is suppressed 
by decreasing temperature T or external magnetic field 
H . The critical behavior of a magnetic impurity shows up 
in the low-T thermodynamics and also low-H properties 
|l7|-po|. Here we would like to determine the latter for 
T = 0. For that purpose, we calculate the spin suscepti- 
bility xo of the impurity at H = and the magnetization 
of the impurity for various values of H. The former gives 
rise to the energy scale (the Kondo temperature) Tk or 
alternatively the screening length £k (xo 1/2V oc £,k)- 
The field-dependent magnetization M(H) exhibits the 
critical behavior of the system, since it contains a term 
proportional to (H/Tk)" with 2 < a < 3 for small H . 
This can be used to distinguish between a local Fermi 
liquid @Jlf[l or an anomalous response case ||. As the 
field is applied exclusively to the impurity spin, we add 
H m = —HS Z to Hq. By using the Hcllmann-Feynman 
theorem, we have M(H) = ($ (H)\ S Z \$ (H)) where 
\^>o(H)) is the ground state in the presence of Hm, and 
Xo = [9M(H)/dH] H _ . xo is numerically evaluated as 
the coefficient of the linear term in an 7i-expansion of 
M(H). For given J and U, the values of H is chosen so 
small, for instance H = 0.00015 at J = 0.5 and U = 20, 
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that the contributions from higher orders in H are neg- 
ligible in comparison with numerical errors (see below). 

The application of the DMRG method to inhomoge- 
neous systems as considered here is not straightforward 
because local couplings must be properly renormalizcd 
by DMRG procedures pl[ . Let us discuss briefly some 
essential points involved in our studies. The number of 
states kept mostly is between 256 and 512. The trunca- 
tion errors are the order of 1CP 10 . In order to work with 
a non-degenerate ground state, our system consists of 
one impurity and an odd number of conduction electrons 
N. The filling factor v = N/L (L: the number of sites) 
equals one, for which DMRG calculations are performed 
with odd L under open boundary conditions, if not stated 
otherwise, but values of v = 1/2 and 3/4 and correspond- 
ingly an even number of sites L are also considered with 
periodic boundary conditions. When L = 4k + 1 (k be- 
ing an integer) and the impurity is located to the centre 
of the chain the ground state has total spin S tot = 
p2j. In this case, the two sites which are added in each 
DMRG step have maximum distance to the impurity, 
and at least three sweeps are taken when the calcula- 
tions are performed. In the presence of Hm , the ground 
state is again taken in the 5'°' = subspace. System- 
atic relative-errors based on truncations are estimated to 
be the order of 10~ 7 . DMRG calculations are done for 
chain lengths 33 < L < 105 so that independent extrapo- 
lations show relative deviations between 10~ 4 and 10~ 3 . 
We consider those as the systematic relative-errors of our 
final results. 

The spin of a Kondo ion is compensated by a spin 
screening cloud surrounding the ion fl23fl . When U = 0, 
the ID system is metallic. The screening length is 
then related to the Kondo temperature Tk via ~ 
vf/Tk where vf is the Fermi velocity. When Jp <C W 
where W is the band width and p = l/2w is the density 
of state at the Fermi level (t = 1 in our calculation) , it is 
£k w ae 1 / Jp /y/Jp. For Tk ~ 10K, the screening length 
extends over thousand lattice spacings a. The strong cor- 
relations between conduction electrons decrease £k sub- 
stantially. Figs. 1(a), and (b) show xo ~ £,k(J) [also 
Xq 1 ~ Tk(J)] for C/ = 0, and 20, respectively, i) a 
crossover: Consider first Fig. 1(a). In order to reach the 
thermodynamic limit, we must have L 3> £,k/cl. This 
limits us to J > 2. The figure suggests a crossover from 
a strong coupling to a weak coupling at J « 2.5. In 
Fig. 1(b), the same crossover takes place at approxi- 
mately J ~ 0.4. ii) a weak coupling regime: Note that 
for J = 0.3 we find xo = 3.73 for U = 20 while for U = 
we obtain xo ~ 10 5 p3| ]. iii) a strong coupling regime: 
One also notices from Fig. 1(a), and (b) that for large 
values of J it is Tk ~ J. The linear dependence sets in 
for U = 20 at much lower values of J (J ~ 1) than for 
U = (J ~ 40). These features on the energy scale Tk 
support previous findings J^-M&0|. 

In Fig. 2, we shows xo as a function of U for J = 4. 



For U >• W, the exchange coupling between sites in the 
chain J ex = 1/4J7 becomes much smaller than J, imply- 
ing that a singlet which extends over a few lattice spac- 
ings is formed between the spins of the magnetic ion and 
the chain, xo saturate at 1/2 J as U ^> W, which cor- 
responds to a local singlet between the spins of the ion 
and the lattice site 0. The situation differs in the case of 
electron attraction (U < 0). One notices a maximum in 
Xo at U c s» -1.0. For U c ( J) < U < the effective Kondo 
coupling is reduced by the attractive electron-electron 
interaction, and £k is enhanced. With increasing attrac- 
tion electrons form more and more on-site pairs. When 
\U\ becomes much larger than the binding energy of the 
Kondo singlet and than W, all the electrons are paired 
except one (remember that N is an odd number). This 
unpaired electron forms a singlet with the magnetic im- 
purity so that xo saturate in the limit U — > — oo [ p4[ . 
Moreover, close to U = 0, the susceptibility xo is linear 
in U, which confirms perturbation results for Tk [ p5| . 
When both Jp and \U\ are much smaller than W, one 
can approximately transform the problem into one of an 
impurity coupled to free electron with an effective Kondo 
coupling given by J**f = J + U/i and . = J. It is i) 
J z eff > J^ f} > when U > Q; ii) < J z eff < J^ {f 
when U c < U < 0; iii) J* ff < and J e i f/ > when 
U < U c . This classification explains qualitatively the be- 
havior of xo in Fig- 2 in the relevant range of U. To 
explore the behavior of xo around U c , we have calculated 
several values of U near U c . As shown in Fig. 2, the cur- 
vature close to U c reveals, with a relative error ~ 10~ 4 , 
that U c is a crossover rather than a critical point. 

It is elucidating to study in more details the different 
states between the impurity and the site 0, which appear 
in the | v I'o(^ = 0)) as well as their weights. They are con- 
structed from the impurity-spin states {|T), II)} and the 
spin states { 1 0) , |t) ||) , |tl) } of site 0, from which we 
can form a singlet \ij} s ) = ^(|T) II) - ll)o lt>); a triplet 
= {|t>olt>, ^(IT>oll> + U>olt», IDoll)}; ^d a 
quadruplet \^ q ) = {|TI) |T), |U) II), |0>„ |T), |0) ||)}. 
We extract the these states from |^o(^ = 0)) by mak- 
ing use of a reduced density matrix. The corresponding 
probabilities are P s , Pt (for each of the three states) and 
P q . When U < U c , one finds P s > P t and P q > P u 
and also has P s ^> P q if J is a sufficiently large. The fol- 
lowing results are found when U > 0: i) for U ~ W, one 
finds P s > P q > Pu h) for U~>W,it is P s > P t > P q . 
When in addition, J ex ^> J, we obtain P s ~ P t ^> P q , 
indicating that the impurity spin is almost fully polar- 
ized by the strong correlated electrons in the chain. In 
contrast, when J ex <C J, we find P s 3> Pt S> P q , imply- 
ing that a local Kondo singlet is formed. Two examples 
are given in Fig. 3 for U = 20 and J = 0.0005, 0.5 where 
also correlation functions (^o(H = 0)| S^sf |*o(-^ = 0)) 
are shown as a function of i. The correlations are long 
ranged, particularly for small values of J. 
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The magnetization of the impurity is well suited for 
studying the critical behavior of the system under inves- 
tigation. Shown in Fig. 4 is the impurity magnetization 
M as functional XoH (~ H/Tk) for various values of 
J and U. The solid line corresponds to the expected 
behavior in the strong coupling limit, i.e., M(\qH) = 
XoH/y/l + A(xqH) 2 . One notices that the data fall onto 
a universal curve with slight but systematic deviations for 
XoH ~ 10 _1 (see inset). Two different regimes are clearly 
distinguishable, a strong coupling regime for lower values 
of XoH and a weak coupling regime in which the magne- 
tization is saturated at M = 1/2. [[l7|-p"9[ In the strong 
coupling regime the impurity spin is well compensated by 
the spins of the electrons in the chain, while in the weak 
coupling regime the singlet is broken by H > Xq 1 - The 
deviations are natural even for {7 = 0, since the calcu- 
lations are performed in real space without any assump- 
tion on the density of states and the values of J are not 
much smaller than W |2(| . This however does not change 
the universality class. Particularly, for small values of 
XoH < 10 _1 , i.e., the critical behavior, M(xoH) behaves 
in the same way for positive and negative U and small 
and large J values at v = 1. For the case of v ^ 1, let 
us consider the lower field behavior quantitatively. Ac- 
cording to ||, one might expect that for small fields H, 
M(xoH) = XoH + ai(xoH) 1 / Kp+1 + higher orders in H. 
When U is finite and v ^ 1, one has 1/2 < K p < 1. 
It turns out that x"(xoH) = d 3 M/d(x H) 3 should be- 
come singular at H = if the above conjecture as re- 
gards M(xoH) holds. In this case, a scaling analysis 
is valid. For finite L, the original question becomes 
whether x'o(Xo(L)H) oc [x {L)H] a( - K "^ with or with- 
out the size-dependent K p (L) for the Hubbard model, 
i.e., Xo (xo(L)H) is divergent as xo{L)H — > if and only 
if K p {L) appear in a(K p (L)). We have computed x"(0) 
for v — 1/2, 3/4, and 1 by the exact diagonalization 
for the length L up to 12. x"(XoH) is accurately evalu- 
ated order by order in a numerical way with the use of 
~ (f(x) — f(x — Sx))/5x for sufficiently small x — 
Xo{L)H and 5x = xo{L)8H. For instance, at v = 3/4 
(L = 12), U = 0.5, and J = 0.01, X "(x) = -12.000173, 
-12.000157, -11.999287 -11.996831, and -11.971354 
for x = 1.476262 x (10" 5 , 10~ 4 , 10~ 3 , 2 x 10~ 3 , 6 x 10~ 3 ), 
respectively, and Sx — 6.6020 x 10~ 5 . We found that 
x"(xoH) = —12.00 at H = is independent of the val- 
ues of v, U (U = 0, 0.5, 10) and J (J = 0.01, 0.5, 
5). For a larger system of L = 33 and N = 25, we 
obtained x"(XoH = 0) = -11.5 by the DMRG calcu- 
lation with keeping 800 states and nine sweeps. Note 
that although this result is less accurate than that given 
by the exact diagonlaization, it is still finite. The above 
analysis turns out that a(K p (L)) equals to zero exactly 
in the same way as for the case U = [pTj. Therefore, 
M(xoH) = XoH - 2(xoH) 3 + higher orders in H as de- 
scribed by a local Fermi liquid. 



In conclusion, we have studied the ground state prop- 
erties of a magnetic impurity coupled to an interacting 
ID system of conduction electrons. We found that a local 
Fermi liquid picture is still valid but that the character- 
istic energy scale, i.e., Tk is substantially enhanced by 
the strong repulsive interaction in the chain and affected 
non-monitonically in the case of electron attraction. We 
infer the validity of a local Fermi liquid description from 
the fact that the magnetization M(xoH) shows the same 
critical behavior for U — and for U ^ 0. 
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FIG. 3. Correlation functions between the impurity spin 
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